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Broken-symmetry quantum Hall (QH) states with filling factors v — 0, ±1, ±2, ±3 in the lowest 
Landau level in bilayer graphene are analyzed by solving the gap equation in the random phase 
approximation. It is shown that in the plane of electric and magnetic fields, the critical line, which 
separates the spin and layer polarized phases at v = 0, extends to the v — ±1 QH states. The 
amplitudes of the gaps in the v = ±1, ±3, and v = ±2 QH states are significantly smaller than the 
C~| ■ amplitude of the v — gap, due to the separate filling of the n — and n = 1 orbital Landau levels 

and the negative contribution of the Hartree term, respectively. It is shown that those values of the 
external electric field where the conductance is not quantized correspond to the minima of the gaps. 

G\ ■ 

■ PACS numbers: 81.05.ue, 73.43.-f, 73.43.Cd 

r rjj 1 I. INTRODUCTION 

U ■ 

' Bilayer graphene is a new material with unique properties.— ~— The possibility of inducing a tunable bandgap by 
. top-bottom gates voltage makes it very promising for applications in various electronic devices. 

Recent experiments in bilayer graphene^— revealed the generation of energy gaps in a magnetic field with complete 
lifting of the eightfold degeneracy in the zero energy (lowest) Landau level (LLL), which leads to new quantum Hall 
states with filling factors v = 0,±1,±2, ±3. While in Refs. HOHEl suspended bilayer graphene was used, bilayer 
graphene samples deposited on SiG^/Si substrates were utilized in Refs. 0H. Because suspended bilayer graphene 
is much cleaner than that on a substrate, the new quantum Hall states in the former start to develop at essentially 
| smaller magnetic fields than in the latter. The theory of the quantum Hall (QH) effect in bilayer graphene has been 

1 1 ! studied in Refs. [U® 

' It was revealed in Ref.[^that the energy gaps scale linearly with magnetic field B in bilayer graphene. It is unlike the 
Kj,. , case of monolayer graphene where a \[~B scaling for the gaps takes place d 1 ' 22 As was suggested in Refs. a strong 

screening produced by the Coulomb interaction is responsible for this modification of the scaling in bilayer; 2 ^ The 
\ physics underlying this effect is the following] 16 ' 17 Due to a nonrelativistic like dispersion relations for quasiparticlesji 
• the polarization function in bilayer graphene in a magnetic field is strongly enhanced as compared to the case of 
| monolayer graphene. In particular, it is proportional to the large mass of the quasiparticles m ~ 10 _2 ™ e . Such a 
qq ■ strong screening radically changes the form of the interaction and leads to the linear scaling. 

Another interesting phenomenon in the v = QH state in bilayer graphene is the phase transition between the 
. spin polarized (ferromagnetic) phase and the layer polarized one in the -B-E_L-plane, where Ej_ is an electric field 
orthogonal to the bilayer planes. It was analyzed in theoretical studies in Refs. I16H19I and observed in experiments in 

> : Refs. tans 

k>( ' ^-ke mam conclusion of the analysis in Refs. [I6IIT7I was tha t these two phases are separated by the critical line 
r> ' Ao = jiB B + e 2 d/l 2 , where [1b is the Bohr magneton, I = y/hc/eB is the magnetic length, Ao = eE±d/2 is the 
jf^ ■ top-bottom gates voltage imbalance, and d ~ 0.35nm is the distance between the graphene layers. This critical line 
is qualitatively consistent with the experimental results in Refs. MIMTol. 

Only the v = QH state was analyzed in Refs. Il6l - fl9l As to the QH states with v = ±1,±2,±3, the main 
experimental results obtained in Refs. 0IH are the following: a) For v = ±1, there are two phases separated by 
approximately the same critical line as that in the v = QH state, b) There is only one phase for v = ±2 and v = ±3 
QH states, c) The v = gap is approximately 30 — 40% larger than the v = ±2 one and significantly (by factor 
10) exceeds the v = ±1 gap. d) For v = 0, v = ±1, v = ±2, and v = ±3, the conductance is quantized except at 
particular values of the electric field E±_ . 

In this paper we extend the analysis of Refs. [l6lfL7l beyond the neutral point and describe the QH states with 
the filling factors v = ±1,±2,±3. As will be shown below, the results of the analysis reproduce correctly the main 
characteristics of the experimental data. 

The paper is organized as follows. In Sec. UH the Hamiltonian of the model, its symmetries, and order parameters 
are described. In Sec. IIII1 by using the Baym-Kadanoff formalism, the gap equation for the quasiparticle propagator 
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including the polarization function is derived. In Sec. IIV1 the properties of the solutions of the gap equation and the 
phase diagram of the LLL are described. In Sec. [Vl we compare our results with experiment. In Sec. I VII the main 
results of the paper are summarized. In Appendix we comment on the bare splitting of the Landau levels with orbital 
indices and 1 induced by a bias electric field. 

II. MODEL 

We will utilize the same model for describing the low-energy electronic excitations as in Refs. 

It can be 

derived from the more accurate four-band model, which takes into account the sombrero shape of the band structure 
for zero magnetic field.— The free part of the effective low energy Hamiltonian of bilaycr graphene is>^ 

#o = ~f d 2 x* Vs (x) ( £ {nt J 2 ) *y.(x), (1) 
where n = p Xl + ip X2 and the canonical momentum p = — ifi,V + eA/c includes the vector potential A corresponding 

2 

to the external magnetic field B. Without magnetic field, this Hamiltonian generates the spectrum E = ±#— , where 
m = ~fi/2vp with the Fermi velocity vp ~ c/300 and 71 w 0.34 — 0.40 eV. The two component spinor field *f>y s carries 
the valley (V = K,K') and spin (s = +, — ) indices. We will use the standard convention: ^J :s = (iPai^B2)k s 
whereas ^k's = {'4 } B2 1 '^ai)k' s- Here A\ and B 2 correspond to those sublattices in the layers 1 (top) and 2 (bottom), 
respectively, which, according to Bcrnal (A 2 — B{) stacking, are relevant for the low energy dynamics. 

The Zeeman and Coulomb interactions plus a top-bottom gates voltage imbalance Ao in bilayer graphene are 
described as (henceforth we will omit indices V and s in the field $y s ): 

ffint = zjd 2 x^ + (x)a 3 ^(x) + 1 Jd 2 xd 2 x' {V(x - x') [ Pl {x) Pl {x') + p 2 (x)p 2 (x')}+ 2V 12 (x - x , )p 1 (x)p 2 (x')} 

+ A J d 2 x<S> + (x)^(x) ~^-J d 2 x*+{x)t ( ^ ) *(a0 , (2) 

where tr 3 is a spin matrix, Z = psB = 0.67 B[T]K is the Zeeman energy, the Pauli matrix r 3 in the voltage imbalance 
term acts on layer components, and £ = ±1 for the valleys K and K' , respectively. The potential V{x) describes 
the intralayer interactions and, therefore, coincides with the bare potential in monolayer graphene whose Fourier 
transform is given by V(k) = 2we 2 /Kk, where k is the dielectric constant. The potential V12 describes the interlayer 
electron interactions. Its Fourier transform is Vi 2 (k) = (2ire 2 / n)(e~ kd / k) . The two-dimensional charge densities P\(x) 
and p 2 (x) are: 

pi(x) = * + (x)P 1 *(x) , p 2 (x) = y+(x)P 2 y{x) , (3) 

where Pi = (1 + £t 3 )/2 and P 2 = (1 — £t 3 )/2 are projectors on states in the layers 1 and 2, respectively. When 
the polarization effects are taken into account, the potentials V(x) and Vi 2 (x) are replaced by effective interactions 
V c g(x) and V\ 2c s{x). The last term in Eq.Q leads to the splitting of the Landau levels with orbital indices n = and 
n = l i 1 ! 4 ! 24 As is discussed in Appendix below, for the values of the magnetic field B < 10T, which are relevant for 
the experiment with suspended bilayer graphene^ and the present analysis, the contribution of this term is small. 
On the other hand, it could be relevant for stronger magnetic fields. Because of that, this term is omitted in the 
analysis of the gap equations in the main text, although we discuss some interesting features this term could lead to 
in the case of stronger magnetic fields in Sees. IIV Al - ITV Dl and Appendix. 

If both the Zeeman and Ao terms are ignored, the Hamiltonian H = Ho+Hi nt , with Ho and H lnt in Eqs. ([T]) and (J2J, 
possesses the symmetry G = U^ K ' 1 (2)5 x U^ K '\2)s^Z^y x Z^y, where U^ v \2)s defines the U{2) spin transformations 
in a fixed valley V = K,K', and Z^y describes the valley transformation £ — >• — £ for a fixed spin s = ± . 16 ' 17 The 
Zeeman interaction lowers this symmetry down to G 2 = U { - K) {1) + x U^ K \l)^ x U ( - K ">(1) + x [/< K ')(1)_ x Z$ x Z^y\ 
where U^ v \l) s is the U(l) transformation for fixed values of both valley and spin. Including the Ao term lowers the 
G 2 symmetry further down to the G 2 = U^ K \l)+ x U^(l)_ x U i * K '\l)+ x U^ K '\l)^. 

The dynamics in the integer QH effect in bilaycr graphene is intimately connected with dynamical breakdown of the 
G and G 2 symmetries. Two sets of the order parameters describing their breakdown were considered in Refs. I1IIT7I . 
The first set consists of the quantum Hall ferromagnetism (QHF) order parameters^ 

p s : = (lpK AlS lpKA lS + ^K, AlS ^K'A lS + ll>K B23 ll>KB 2S + B^K- B % s) , (4) 

jl s : = (lpK AlS lpKA lS + ^ KB2S ^KB 2 s ~ IpK'B^K 'B 2 s ~ ^K'A^K >A 1S ) ■ (5) 
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The order parameter ([4| is the charge density for a fixed spin whereas the order parameter ([5]) determines the 
charge-density imbalance between the two valleys. The corresponding chemical potentials are /i s and p, s , respectively. 

The second set consists of the magnetic catalysis (MC) order parameters,— i.e., the Dirac A s and Haldane A s mass 
terms: 

A s : (^ f s T 3 ^ s ) = 0'KA 1 s^KA 1 s-tpKB 2S i J KB 2S + ^K'B 2S i'K'B 2 s-tpK l A lS ' l ( J K'A 1 s) , (6) 
As : (*^T 3 * S ) = {^kA 1 s^KA 1 s+^k, Ais ^K'A 1 s-^kb 2S ^KB 2 s~'^k, B2S ^K'B 2 s) ■ (7) 

Clearly, the order parameter (j6|) describes a charge density wave in both the K and K' valleys. While this order 
parameter preserves the G 2 symmetry, it is odd under time reversal.— On the other hand, the order parameter (|7|). 
connected with the conventional Dirac mass A s , determines the charge-density imbalance between the two layers X 
Like the QHF order parameter ([5]), this mass term completely breaks the Z^y symmetry and is even under time 
reversal. It is important that in both monolayer and bilayer graphene, these two sets of the order parameters 
necessarily coexis t 16 ' 17 ' 28 and are produced even at the weakest repulsive interactions between electrons (magnetic 
catalysis 2 ^—). The essence of this phenomenon is an effective reduction by two units of the spatial dimension in the 
electron-hole pairing in the LLL with energy E = . 

Let us also emphasize that unlike a spontaneous breakdown of continuous symmetries, a spontaneous breakdown 
of the discrete valley symmetry Z^y, with the order parameters ( , 3 / ^ v E' s ) and ( , I / J£t3 , I' s ), is not forbidden by the 
Mermin- Wagner theorem at finite temperatures in a planar system.— Also, because the valley and layer indices are 
equivalent in the LLLji these order parameters also describe a breakdown of the symmetry between the top and 
bottom layers. 

Note that because of the Zeeman interaction, the SU^ v '(2)s is explicitly broken, leading to a spin gap. This gap 
could be dynamically strongly enhanced^. In that case, a quasispontaneous breakdown of the SU^ v \2)s takes place. 
The corresponding ferromagnetic phase is described by the chemical potential ^3 = (/Lt+ — /i_)/2, corresponding to 
the QHF order parameter (^a^), and by the mass A3 = (A+ — A_)/2 corresponding to the MC order parameter 

Recall also that in bilayer graphene, the LLL includes both the n = and n = 1 LLs, if the Coulomb interaction 
is ignored^ Therefore, in the LLL approximation, there is an approximate orbital symmetry Z2L connected with a 
n = — > n = 1 transformation. 



III. GAP EQUATION 

The gap equation for the full propagator G in the LLL approximation in bilayer graphene was derived in Ref. ITtI 
In the derivation, the Baym-Kadanoff (BK) formalism^ was utilized. In this section, we describe the main features of 
the gap equation. Note that because a characteristic scale in the bilayer dynamics in a magnetic field is the cyclotron 
energy huj c ~ 2.19P[T]meV, the applicability of the LLL approximation implies that the LLL energy gaps should be 
smaller than Huj c . As we will sec, this condition is fulfilled in bilayer graphene. 

We will analyze the gap equation for the full quasiparticle propagator G with the order parameters introduced in Eqs. 
(U])-([7]). The mean field approximation with the polarization function calculated in the random phase approximation 
(RPA) will be used. The corresponding two-loop BK effective action is a functional of the full quasiparticle propagator 
G and it has the formi^I 

T(G) =r n (G)+F 2 (G), T (G) = -iTi-[LnG- 1 +S- 1 G-l] , (8) 

where S is the free Green's function, Tq is the free part, and r 2 takes into account the interaction effects in the first 
order of perturbation theory, 

T 2 (G) = - J d 3 ud 3 u' |^tr[G( it, u')G(u', u)] V e f[(u — u') + tr [Pi G(u, u') P 2 G(u' , u) ] Vil(u — u') 

- itr [G(u, u)} tr [G(u', u')] V cS (u - «') - tr [ Pi G(u, u)} tr [ P 2 G(u', u') ] V IL (u - u')| ■ (9) 

Here u = (t,r),t is the time coordinate, r = (x, y), and Vn,(u) = Vi2eff(w)— V e s(u) is a combination of the interlayer and 
intralayer interactions (recall that the polarization contributions are included in the potentials V e s(u) and Vi2 e ff(u)). 
Note that while here the trace Tr, the logarithm, and the product 5 _1 G are taken in the functional sense, the trace 
tr runs over layer, valley and spin indices. The stationary condition ST/8G = leads to the gap equation for the 
quasiparticle propagator G in mean field approximation. 
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We will use the Landau gauge for a two dimensional vector potential, A = (0, Bx), where B is the component of 
the magnetic field B orthogonal to the xy plane of graphene. Then, the free Green's function S(ui, u 2 ) can be written 
as a product of a translation invariant part S(ui — u 2 ) times the Schwinger phase factor— 



S(ui, u 2 ) = cxp 



. (xi + £2X2/1 



2l 2 



S(ui 



U2) 



(10) 



Similarly, we can separate the Schwinger phase factor from the translation invariant part in the full propagator 

\xi +£2X2/1 -2/2)' 



G(ui,u 2 ) = exp 



2P 



G{ui - u 2 ) . 



(11) 



For the LLL with the orbital numbers n ~ 0, 1, the translation invariant part of the time Fourier transform of the 
free propagator takes a simple form: 



1 



with 



2ttZ 2 



{ 2P 



1 



uj + /.io — sZ + £Aq + iSsgnuj 



(12) 



(13) 



where /io is the electron chemical potential, Z is the Zccman energy, and the projector P_ = (1 — Tz)/2. 

Motivated by expression (|12p for S^ s (r; uj), we will use the following ansatz for the full propagator with the param- 
eters n s {n), p. s {n), A s (n), and A s (n) related to the order parameters in Eqs. (JU) - (JZ]): 



G{ s (r;w) 



2nl' 



■ cxp 



r 

472 



r 

272 



G f i s (w)ii 



r 

2P 



(14) 



where 



UJ — E^ na + iSsgTLUJ 



and 



E Sns = -(Vs(n) + A s (n)) + £(fi s (ri) - A s (n)), n = 0, 1, 



(15) 



(16) 



arc the energies of the LLL states depending on the parameters /i s (n), /i s (n), A s (n), andA s (n). Note that because for 
the LLL states only the component tpB 2 s (^Pa 1 s) of the wave function at the K(K') valley is nonzero, their energies 
depend only on the eight independent combinations of the QHF and MC parameters shown in Eq. (fTB"]) ^ Our goal 
is to find the energies E^ ns from the gap equation. 

It is convenient to solve the gap equation in momentum space. Then, the Fourier transform V e ff(uj, k) of V c s{u) 



is 
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V cS {uj, k) 



2ne z 



1 



fc + ^fn( w ,k 2 ) 



(17) 



with II(w,fc 2 ) = nii(aj,k.) + IIi2(w,k.), where the polarization function Lly describes electron densities correlations 
on the layers i and j in a magnetic field. For the polarization function n(w, k 2 ) we use the expression in the RPA 
approximation modified by the presence of a quasiparticle bare gap A term (see Eq.(A24) in Ref. flTi ). As to the 
potential Vil(u), the gap equation contains its Fourier transform only at zero frequency and momentum^ 



V Ih (cj = 0,k = 0) 



2ne 2 d 



(18) 



The inclusion of the polarization effects is crucial for ensuring the linear scaling in bilayer graphene. We utilize 
the frequency independent order parameters /i,/2, A, A and the static approximation for the polarization function 
will be used, n(a;,k 2 ) — > 11(0, k 2 )^ The explicit expression for the polarization function calculated in the RPA 
approximation can be found in Ref. [TjJ. 
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Using the ansatz (fl~4)) . one finds that the gap equation 5T/SG = leads to the following system of equations for 
the energies E^ sn (for details, see Ref. [l7j): 



Ho - sZ + £A 
1 



[sgn (E (0s ) h{x) + sgn (E^ ls ) I 2 {x)} 



Airl 2 



(A 1 +A 2 )V cS (0) + 



^A 2 + ^A 1 1 V 1L (0) 



Ho - sZ + £A - 
1 



2ml 2 

(A 1 +A 2 )V eS (0) 



2 2 
gn (E i0s ) I 2 {x) + sgn (E ils ) I 3 (x)] 



A 2 



V5l(0) 



where the integrals h(x) are 



dyfi(y)e v 
o K^/^y + 47rn(y) 



(19) 



(20) 



(21) 



with fi(y) = (1, y, (1 — y) ) for i = 1,2,3, respectively. Here the dimcnsionless variable x = 2h /e m I = 
(4hoj c /a 2 ji)(vF/c) 2 ~ 0.003B[T], where a = 1/137 is the fine-structure constant and we used the values 71 = 0.39eV, 
hw c = h 2 /ml 2 = 2.19B[T]meV, v F = 8.0 x 10 5 m/s (see Ref.[l|). The quantities A x and A 2 are Ai = £ sgn( E_ lns ) , 
^2 = E„, s sgn(^i„ s ). 

The terms in the first and second square brackets on the right hand side of Eqs. (fT9"| and (|2"U)) describe the exchange 
and Hartree contributions, respectively In the latter, the terms with the factor (A\ + A 2 )/4irl 2 are proportional to 
tr[G(0)]ji£ i.e., the density of charge carriers. There are other Hartree contributions, such as those taking into account 
the background charge of ions in graphene, the charge in the substrate, etc., which are not included in the equations. 
Due to the overall neutrality of the system, all these contributions should exactly cancel. As a result, only the Hartree 
terms with the factor (A\ — A 2 )/4irl 2 survive. They describe the capacitor like interactions in bilayer graphene. Thus, 
the final equations, taking into account the neutrality condition, are: 



E^os = -Ho + sZ - £A 

E(i s = -/-to + sZ - £A 
The filling factor 



[sgn (E (0s ) I x (ar) + sgn (E ils ) I 2 (x)} 



2ml 2 

h 2 

^—^ [sgn (E (0s ) I 2 {x) + sgn (E (u ) h(x)] 



" = -^sgn(S e „ s ) = --{A 1 +A 2 ) 

tylS 



8irl 2 

J_ 

8nl 



(A 1 - A 2 )V lh (Q) 
2 (A 1 -A 2 )V Uj (0) 



(22) 
(23) 

(24) 



takes values 0, ±1, ±2, ±3±4. In general, there are many solutions of the gap equations ([22]) and (|23|) at a fixed filling 
factor v. The density of the thermodynamical potential of the system for each solution can be calculated by using 
the BK effective action © and the fact that solutions are extrema of this action. We find 



n = -i 



8tt 2 / 2 



E 



w — Mo + sZ — £A 



- 1 



8ttZ 



lo — E£ ns + iSsgnuj 
|pEEE [ E t™ + sZ £Ao 



£=± s=± n=0,l 



sgn(£f„ s ) . 



(25) 



IV. SOLUTIONS AND PHASE DIAGRAM IN THE LOWEST LANDAU LEVEL 



The chemical potential is one of thermodynamic variables of the thermodynamic potential (|25|) . However, in 
experiments with bilayer graphene, the top and/or bottom gates control the charge density (filling factor v) rather 
than the chemical potential. Therefore, it is convenient to perform a Legendre transform and use the charge density 
as a thermodynamic variable. 

Usually, there is a one-to-one correspondence between the chemical potential of a system and its charge density. 
However, in the case under consideration, such a correspondence is absent. Indeed, since we consider an ideal system 
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at zero temperature and zero LLL width, the chemical potential can take arbitrary values within the gap between the 
filled and unfilled levels. Since the LLL states corresponding to the filling factors v = 0,±1,±2, and ±3 are gapped, 
this is a general situation in the present dynamics. 

Actually, this simplifies the analysis of these states: because as soon as the filling factor is fixed, the chemical 
potential becomes an irrelevant parameter and does not affect the energy density of the system for any solution 
studied. Explicitly, using the charge density 



P 



/-\-oo J -i 
— tr[G( W ;0)] = -— E E E sgn(i^„ s ) (26) 
-oo ~* r l r. L — n i 



£=± s=± n=0,l 

and performing the Legendre transform, we find the following free energy density: 

s = n + ^op = -g^j E E E [ E t»> + 1 10 + sZ - s s n ( E s™) ■ ( 27 ) 

f=± s=± n=0,l 

Since — /iq additively enters E^ ns defined in Eqs. ([22|) and (|23l) . it is clear that /iq cancels out in ([27| for all solutions 
with a given filling factor. 

Therefore, our strategy in the analysis of the gap equations is the following one: At each fixed filling factor v and 
given values of the controlling parameters B and E±, we find all possible solutions of Eqs. (|2"2"|) and (|2"3"]l with a 
chemical potential allowed by the filling v, and then determine the ground state as the solution with the lowest free 
energy density (|2"T|) . 

The gap equations (|22j) and ([23| form a system of algebraic equations for the eight energies E^os and E^\ s of 
quasiparticle levels. The solutions having the same filling factor differ in the signs of particular energies but have the 
same combined values of A\ + A2 (and consequently v). Therefore in order to describe the solutions corresponding 
to different filling factors v, it is sufficiently to specify only the signs of their energies. 

By pre-setting the signs of the energies, one finds that for v = 0, ±1, ±2, ±3, and ±4, there are (®) = 70, 
(3) = 56, (2) = 28, (^) = 8, and (T) = 1 solutions, respectively. Therefore there are 256 solutions all together. In 
what follows, we will consider the solutions only with filling factors from to +4 (the solutions with negative values of 
v can be obtained from these ones by choosing opposite values of the signs of the energies) . This reduces the number 
of the solutions to 163. For concreteness, we also choose positive values for a magnetic field, B > 0, and nonnegative 
values for an applied electric field, E±_ > 0. 

As is discussed below, the analysis of Eqs. (|22|) . ([23)) and (|27|) with a chemical potential allowed by the filling iS^L 
shows that for each filling factor, there are only one or two major solutions that have the lowest energy density and 
therefore are relevant for the phase diagram. Only for v = 1 and v = 3, and in a very narrow region of the values of 
B and E± (where the major phases do not exist), there might exist some additional (marginal) phases. 

In this section, the properties of the solutions of gap equations (|2"2"|) - ([2"3"|) will be described. The comparison of these 
solutions with experiment will be considered in Sec. [V] 

A. The v = QH state 

The solutions of the gap equations at the neutral point were analyzed in Refs. fl6lli~7l ]. Here we briefly describe 
those results. 

There are two competing solutions: I) the ferromagnetic (or spin polarized (SP)) solution and II) the ferroelectric 
(or layer polarized (LP)) solution. The signs of the energies E^ ns in the ferromagnetic solution are negative for the 
states with a spin opposite to magnetic field and positive for the states with a spin along magnetic field: 

sgn(Si, ,-) = -1, sgn(£Li i0 ,-) = -1, sgn(£i,i _) = -1, sgn (£U,i.-) = -1, 

sgn (£1,0,+) = 1, sgn (£_i,o,+) = 1, sgn (E 1A ,+) = 1, sgn (£Li,i,+) = 1 . (28) 

In the LP solution, the signs of the energies are correlated with the signs of the valley index £ (recall that the valley 
and layer indices are equivalent in the LLL): 

sgn (£1,0,-) = -1, sgn(£a, , + ) = -1, sgn (E 1A _) = -1, sgn(£'i,i j+ ) = -1, 



sgn(£_x,o,-) = 1, sgn (JS_i j0 ,+) = 1, sgn^i^,-) = 1, sgn (£L 1: i i+ ) = 1 



(29) 
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The analytical expressions for the energies E^ ns of these two solutions can be easily found from Eqs. (|22|) and (|23l) 
Then, from Eq. (f2~T)) with fi = 0^1 one finds their energy densities: 



8ml 



:{h+2L 2 + h 



(30) 



and 



£Lp --^l Aa ~7^ ' 8mP 



(h + 2I 2 + h) 



(31) 



where the integrals Ii are defined in Eq. (|2ip and the dielectric constant k is a free parameter in this model. Comparing 



v=0, \i Q =0 [meV], k=2 




20 30 

E x [mV/nm] 



40 



50 



FIG. 1: The phase diagram for the v = QH state in the (E±, B) plane. Here the dielectric constant k = 2. 

the expressions for £sp and £lp, one concludes that there is a critical line separating the SP and LP phases in the 
A -B plane^iil 



A? = Z 



d 2 ' 



(32) 



In order to compare our results with experimental data in Refs. 0-E3, it is convenient to express the top-bottom 
gates voltage imbalance Ao through the electric field, E± = 2Ao/ed. Then, the critical line in the E±B plane takes 
the form 



(33) 



Using Z ~ 0.67B[T]K and I — ^hc/eB, this relation can be rewritten as 



mV 4 4 

E™[— - ~ (0.33 + — B T 
nm k 



(34) 



We plot the phase diagram of the v — state in Fig. [T] (compare with Refs. [l6lfr7l ). The I (II) area is that where 
the SP (LP) solution is favorite. The red dashed line and the B axis compose the boundary of the region where the 
two solutions coexist: the solution I does not exist to the right of the (red) dashed line in the region II. The black 
bold line is a critical line between the phases I and II. For balanced bilayer (E± = 0), the SP solution is favorite due 
to the presence of the Zeeman term. At fixed B for sufficiently large electric field E±, the LP solution is realized. 

The gap A^o for the Hall state v = is equal to the difference between the lowest empty (positive energy) level 
and the highest filled (negative energy) level. For the SP solution, it is 



A„ =0 — Ei 



£-1,1,- = 2 



Z Ei 



2ml'- 



(h + h) 



E ± < —\Z 
ea 
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(35) 
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and 



A„ =0 = E- 



-E 



1,1.+ 



= 2 



h 2 . . 2e 2 d 

— (h+h)--jj- 



2ml 



{h+h) + -^\ (36) 



for the LP solution. 

At a fixed value of the magnetic field, the SP gap decreases with increasing E± 



up to the point E± 



2\Z 



h /2ml (I2 +Is)]/ed, where it vanishes. On the other hand, the LP gap always increases with E± 



the expression in Eq. (j33j) for the critical line, the SP solution is favored for E± < Ef, while at Ej 
solution is realized. At the critical line, the gaps of the two phases coincide and are given by 



/\cr _ 
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(37) 



Let us describe the phase transition between the SP and LP phases in more detail. In Fig. [21 the gap A„ = o and the 
energy spectrum E^ ns are shown as functions of the electric field E± at the fixed value of the magnetic field B = 2 [T] . 
The critical value of E± is Ef ~ 5.04mV/nm and Af =0 ~ 3.73 mcV in this case. 




FIG. 2: The energy spectrum and the gap of the v = QH state as functions of the electric field at the fixed magnetic field 
B — 2 [T]. Here the dielectric constant k — 2. 



As one can see in Fig. [21 ignoring a small splitting of the LLs due to the electric field, the symmetry of the SP 
phase at E± < Ef is G 2 = U^ K \l)+ x U^ K \l)- x U^ K '\l)+ x C/( A ")(1)_ x Z { 2 y } x Z { 2 ~ ] considered in Sec. El Now, 
at the critical value Ef, a jump in the energy spectrum takes place (see Fig. |2j). Since the value of the gap A" =0 
at E± = Ef is nonzero, this leads to a jump transformation of the G2 symmetry into the U^ K '{2)$ x U^ K '(2)s one 
at E± > Ef (the LP phase), if the Zeeman term is ignored. The presence of such a jump suggests that this phase 
transition is a discontinuous (first order) one^ 

A noticeable feature of this transition is a kink singularity in the gap with a minimum at E±_ = Ef clearly 
seen in Fig. [2] As will be discussed in Sec. [V] this fact is important for understanding the behavior of the two- 
terminal conductance in the experiment in Ref. [jj. It is interesting that as shown in Appendix, the last term in the 
Hamiltonian ([2]), which is responsible for the bare splitting of the Landau levels with n — and n = 1, transforms the 
kink singularity into a (stronger) jump one. Although for magnetic fields B < 10T, which are relevant for the present 
analysis, this jump is very small, it could become relevant for dynamics with higher magnetic fields. 

B. The v = 1 QH state 

There are two main solutions at the v = 1 filling factor (see Fig. [3]). The first solution is: 

sgn(£ li0 _) = -1, sgn(E , _i )0 _) = -1, sgn(#i,i_) = -1, sgn (E- lfl ^) = -1, 



sgn(£ M ,+) = -1, sgn (£-1,0,+) = 1, sgn(E hh+ ) = 1, sgn (£-1,1,+) = 1. 



(38) 
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This solution is closely connected with the SP one in Eq. ([28]). It differs from the latter only in that its LL with 
£ = 1, n = 0, s = + is filled. It would be appropriately to call this solution a partially spin polarized (PSP) one. 
The second solution is: 



sgn(£i,o,-) = -1, sgn(£a i0 , + ) = -1, sgn {E 1A _) = -1, sgn(Ei ) i )+ ) = -1 
sgn (E-1,0-) = -1, sgn (£_i )0 ,+) = 1, sgn (£Li, = 1, sgn(S_i,i >+ 



1. 



(39) 



It is closely connected with the LP solution in Eq. (|29|) . The difference is in that the LL with £ = — 1, n = 0, s = 
is now filled. It would be appropriate to call this solution a partially layer polarized (PLP) one. 
The energy densities of the partially filled spin polarized and layer polarized solutions are given by: 



£ 



1 



PSP 



and 



£plp 



2irP 



2nP 



iZ 



2ml 
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zE\d 
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Z 



h 2 lr T . 9e 2 d 3eE ± d 
— {h + h+ h) x 



2ml 



(40) 



(41) 



respectively. Comparing these two energy densities, we conclude that the critical line between the PSP and PLP 
solutions exactly coincides with that between the SP and LP ones (see Eq. (|34]l). It is quite noticeable that this 
result is in complete accord with experimental data in Ref. @ according to which the phase transition between the 
two v = +1 states takes place near the same electric field at which the transition between the two v — states is 
observed. 
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FIG. 3: The energy spectrum and the gap of the v = 1 QH state as functions of the electric field at the fixed magnetic field 
B = 2 [T]. Here the dielectric constant n — 2. 



The gap A„ =1 and the energy spectrum E^ ns for the i> = 1 QH state are shown in Fig. [3] As one can see, a jump in 
the spectrum at E± = E" ~ 5.04mV/nm is similar to that in the v = QH state. The phase transition at E± = E" 
corresponds to the transformation of the U {K) (l) + x U {K \l)- x [/ (K ' } (1) + x [7 (K,) (1)_ x Z { 2 ~ ] symmetry in the PSP 
phase into the U^ K \2)s X U' K '(1)+ x {/^ '(1)- in the PLP one (as in the v = case, we ignore small splittings in 
the spectrum due to the Zeeman term and the electric field). 

As to the energy gap, there is an essential difference in its behavior in comparison with the v = gap. The energy 
gap for the PSP solution is found to be 



A, 



Ei,o. 



2ml' 



in the region E^ ns < E±, where it is the ground state. In the region E$. 
ground state, the energy gap takes the same form: 



2/ 2 ) 

s > E 



_L • 



(42) 



where the PLP solution is the 



E. 
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2ml' 
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As a result, the gap in the ground state is a smooth function of E± (without a kink), unlike the v = case. As will be 
discussed in Sec. El this feature and the phase transition at E± = Ef are important for understanding the behavior 
of the conductance at v = 1 . 

Note that both gap (|4"2"j) and gap (f4"3"f are expressed as differences of the LL energies with different orbital number 
n. This takes place also for the v = 3 QH state (see Sec. IIVDI below). The polarization of the LLs states with 
different orbital numbers at odd-integer filling factors was predicted in ReffTTI based on the Hund's rule. 

The following remark is in order. In addition to the two solutions described above, there are also two other, 
marginal, solutions. The first of them is the same as the PSP solution (|3"5)l except that instead the state with £ = 1, 
n = 0, s = +, it has the state with £ = 1, n = 1, s = + being filled. The second marginal solution is the same as the 
PLP solution (|3T))) except that the state with £ = — 1, n = 1, a = — is filled instead that with £ = — 1, n = 0, s = — . 
These marginal solutions have higher energy density than solutions (|3"8)) and (f3T)]) . however, they exist in a slightly 
larger region. Therefore, the marginal solutions can describe the ground state of the system only in a small region. 
This issue will be considered in detail elsewhere. 

At last, we comment on how the bare splitting of the Landau levels with orbital indices and 1, discussed in 
Appendix, influences the v = 1 solution. First of all, the value of the critical electric field in the v = 1 state remains 
the same as in the v = one. Also, as in the v = case, the bare splitting leads to a jump in the v = 1 gap at the 
critical point. Its value coincides with that in the v = gap. Although the jump is very small for magnetic fields 
B < 10T, it could be relevant for dynamics with stronger magnetic fields. 



C. The v = 2 QH state 



There is only one solution that describes the ground state at the filling factor v = 2 (see Fig. @|: 
sgn(£i, ,-) = -1, sgn(£Li, ,-) = -1, sgn(£a,i _) = -1, sgn(^_i,i_) = -1, 



sgn(£i i0 ,+) = -1, sgn(.E M , 



-1, sgn(£_i i0 , + ) = 1, sgn = 1 



(44) 



It can be obtained from the PSP solution (|35|) in the v = 1 QH state by filling the LLs with £ = 1, n = 1, s = + 
or, alternatively, from the PLP solution (|3T)|) by filling the LLs with £ = — 1, n = 1, s = — . In fact, perhaps more 
clear understanding of this solution can be obtained by considering the SP and LP solutions in the v = QH state. 
For the SP solution (|2"5]l , its LLs are spin polarized and the LLs with s = — are filled. For the LP solution (f2T)]) . its 
LLs are layer polarized and the LLs with £ = 1 arc filled. Solution (j44|) can be obtained from the SP one by filling 
the LLs with £ = 1 or, alternatively, from the LP solution by filling the LL with s = — . It would be appropriate to 
call it a partially spin-layer polarized (PSLP) solution. The energy density of this solution is 

h 2 /T nr T , e 2 d eE ± d~] 

;{I X +2h+h)- 
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The gap A^ = 2 and the energy spectrum E^ ns for the v = 2 QH state are shown in Fig. [4] The gap is 
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at E± < E c l ~ 5.04mV/nm, and 
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(46) 



(47) 



at E ± > E c , 



5|) due to the negative Hartree 
2 energy spectrum. Instead, a 



Note, that at E± = the gap A y= 2 ([41))) is smaller than the gap A„ = o 

contribution. 

As one can see, unlike the v = and v = 1 QH states, there is no jump in the v 
LLs crossing takes place when the electric field E± is between 3.45 mV/nm and 6.41 mV/nm. For all values of E±, 
the symmetry of the v = 2 QH state is U^ K \2)s x U^ K • l (l)+ x U^ K '(1)-, if small splittings in the spectrum due to 
the Zeeman term are ignored. 

As in the v = QH state, there is a kink in the A^ = 2 gap at E± = E" ~ 5.04mV/nm, however, its structure is 
quite different: While there is a minimum in A^o at E± = E™, there is no either minimum or maximum in A^ = 2 at 
this same point (compare Figs. [2]and[3|. As will be discussed in Sec. [V] this fact is important for understanding the 
behavior of the conductance at v = 2. 

The influence of the bare splitting term of the n — and n — 1 LLs on the v — 2 state is described in Appendix. 
This term leads to the crossing value E^° ss that is different from the critical value E" in the v = and v = 1 states, 
although their difference is small for B < 10T. The kink singularity in the v = 2 gap remains unchanged. 
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FIG. 4: The energy spectrum and the gap of the v — 2 QH state as functions of the electric field at the fixed magnetic field 
B = 2 [T]. Here the dielectric constant n — 2. 



D. The !/ = 3 QH state 

There is one main solution and two marginal ones in this case. The latter will be considered elsewhere. As to the 
main solution, it is given by (see Fig. [5]) 



sgn(£i,o,-) = -l,sgn(£i i0 ,+) = -1, sgn (£1,1,-) = -1, sgn (£1,1,+) = -1, 

-l,sgn 0EL M ,_) = -l,sgn(£_ 1)li+ ) = 1 



(48) 
(49) 



sgn (£-1.0,-) = -l,sgn(£ , _i j0)+ ) 

It can be obtained from the PSLP solution (|44|) in the v = 2 QH state by filling the LLs with orbital n = and 
£ — — 1, s — +. Its energy density is 
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The gap A„ = 3 and the energy spectrum E^ ns for the v = 3 QH state are 
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FIG. 5: The energy spectrum and the gap of the v — 3 QH state as functions of the electric field at the fixed magnetic field 
B = 2 [T]. Here the dielectric constant n — 2. 
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(51) 



for all values of E± . 

One can see in Fig. [5] that as in the v = 2 QH state, there is no jump in the v = 3 energy spectrum. Instead, a 
LLs crossing takes place when the electric field E± is between 1.48 mV/nm and 3.45 mV/nm. The symmetry of the 
v = 3 QH state is the same as that of the v = 2 one, U {K) {2) S x £/ (K,) (l)+ x [/( K ')(l)_ (if small splittings in the 
spectrum due to the Zeeman term are ignored). The important difference between these two QH states is that the 
A„ = 3 gap in Eq. (|51[) is expressed as a difference of the LL energies with different orbital number n (similarly to the 
gaps of the v — 1 state). Therefore we will call this solution a partially spin- layer-orbital (PSLOP) polarized one. 

As in the case of the v = 1 QH state, there is no kink in the gap A„=3. This fact is relevant for understanding 
the behavior of the v = 3 conductance discussed in Sec. |VJ As shown in Appendix, the influence of the bare splitting 
term of the n = and n = 1 LLs on the v = 3 state is reduced to adding a small (for B < 10T) term in the gap. 

E. The v = 4 QH state 

The ground state is determined by the following solution in this case: 

sgn(£ , i, ,-) = -1, sgn(£'i,o,+) = -1, sgn(-Ea,i,_) = -1, sgn(£ , M , + ) = -1, 



sgn (E. 



1,0,- 



= -1, sgn(£_i, , + ) = -1, sgn (£_!,!_) = -1, sgn (£-1,1,+) 



-1, 



(52) 



where the last empty LL in the PSLOP solution (|49p (with £ = —1, n = 1, s = +) is now being filled. Thus, the v = 4 
QH state describes the filled LLL. The energy density of this solution equals 



£ v= a — — 



J (7 1 +2/ 2 +7 3 ), 



(53) 



The energy spectrum E^ ns for the v = 4 QH state is shown in Fig. |H1 At = 0, its symmetry is G = U tyK \2)s x 
U( K \2)s x Z^y x Zjy^ (if small splittings in the spectrum due to the Zeeman term are ignored). At nonzero E±, 
the symmetry becomes U^ K '(2)s x U^ K '(2)s- Obviously, this QH state describes the gap between the LLL and the 
n = 2 LL. 



> 



v=4, B=2 [T], n Q =0 [meV], k=2 



-2. 



II 



-10 



»»....SSSSSo sSSS5 
^SgagSK-:--- 




^=+1, n=0, s=- 

^=+1, n=0, s=+ 

- - ^=+1, n =1,s=- 
5=+1, n=1, s=+ 
+ |=-1, n=0, s=- 
« £=-1,n=0, s=+ 
° ^=-1 , n=1 , s=- 
• i;=-1 , n=1 , s=+ 







10 15 

E± [mV/ran] 



20 



25 



FIG. 6: The energy spectrum of the v 
Here the dielectric constant k = 2. 



4 QH state as a function of the electric field at the fixed magnetic field B = 2 [T] 



COMPARISON WITH EXPERIMENT 



In this section, we will compare the results of our analysis with the recent experimental data in Refs. 0-E The 
comparison is hampered by the presence of disorder in real bilayer samples, which is ignored in the present analysis. 
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Still, as will be shown below, the results of the analysis reproduce correctly the main characteristics of the experimental 
data. 

Let us start from the critical line between the SP and LP solutions in Eq. (|34| . One can see that the maximum value 
of its slope is obtained at the smallest permissible value K = 1 for the dielectric constant: it is about 4.7mVnm~ 1 T~ 1 . 
On the other hand, the value of the slope in experiments^^ is about llmVnm -1 T - . This discrepancy may have 
its roots in disorder, depending on a type of the latter. For example, an external electric field is more effective in a 
clean sample (considered in our model) and, therefore, its critical value should be smaller than that in a real sample 
with charged impurities. On the other hand, neutral impurities might act just in opposite direction by diminishing 
the role of the Hartrcc interaction in the gap equations, thus favoring the LP solution. This problem deserves further 
study. 

Let us now turn to the energy gaps as functions of magnetic field at E± = shown in Fig. [7] Comparing this 
figure with Fig. 2 in Ref. H, one can see that the hierarchy of the gaps is the same: A^o > A I/=2 > A„ = i = A„ =3 . 
Moreover, the theoretical and experimental values of A^o are close. The worst description in the model takes place 
for the smallest A v= i gap: while the experimental A„ = i is significantly (by factor 10) less than the experimental 
A v —o, the gap A^i is twice less than A„ = o in our model. This fact can be related to the observation in Ref. [H that 
the QH states with higher filling factors are subject to additional sources of disorder. Then, it is not unreasonable to 
assume that the smallest gap, A u= i, is most sensitive to disorder and therefore suppressed stronger. As to the v = 4 
QH state, it describes a gap between the LLL and the n = 2 LL, which is A^ = 4 = \/2hw C} up to small corrections. In 

E ± =0 [mV/nm], fi =0 [meV], n=2 

20 1 . 1 1 . . 1 




FIG. 7: The energy gaps as functions of magnetic field for different filling factors at zero electric field and k — 2. 

the rest of this section, we compare the behavior of the energy gaps as functions of electric field in Figs. [2] -[5] with 
that of the two-terminal conductance in experiment in Ref. 



A. Energy gaps and conductance: even v 

One of the main experimental results in Ref. is that the two-terminal conductance is quantized except at particular 
values of the electric field E±. Let us start from the states with even filling factors. 

For the v = state, the experimental data are the following. The conductance in the v = state is quantized 
except two values of E±, which increase with B£ Let us argue that these two values correspond to E± = ±E", 
where Ef is the critical value at which the phase transition between the SP and LP phases takes place (see Fig. [2] 
recall that for convenience, we consider only nonncgativc values of E± in the analysis). As one can see in this figure, 
the gap A„ = o has a maximum at E± = and a minimum at E± = Ef. This implies that while the conductivity is 
suppressed both around zero E± and for electric fields larger than Ef, it is enhanced (has a maximum) at E±_ = E". 
This picture agrees with that in Fig. 2C in Ref. 0. 

As to Ef, while its experimen tal val ue is Ef ~ 29mVnur 1 at B = 2.65T,! it is Ef ~ e.emVnm- 1 at B = 2.65T 
and k = 2 in our model (see Sec. IIV A[) . This discrepancy can be traced to the discrepancy between the experimental 
and theoretical values of the slope of the critical line, which in turn could be caused by disorder in real graphene 
samples (see the discussion above). As one can read from Eq. (f34|) . the slope of the critical line in our model is 
2.5mVnm~ 1 T~ 1 at k = 2 while its experimental value is about llmVnm _1 T _1 , i.e., by factor 4.4 larger. Taking this 
factor into account, we find that it will increase Ef from 6.6mVnm _1 to 29mVnm _1 that agrees with the experimental 
value. 
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For the v = 2 state, the picture is quite different. In this case, the conductance quantization is broken only at 
E±_ = OZ It agrees with the following facts (see Fig. a) unlike v = 0, there is only one phase determining the 
ground state for v = 2 ; b) for v = 2 the gap at E± = is approximately 30 — 40% less than that for v = and 
monotonically increases with E± ; c) as already mentioned above, the disorder suppresses the v = 2 gap stronger than 
the v = one. Let us also note that as was observed in Ref. |7], in stronger magnetic fields, B > 7.8T, the v = 2 
conductance is quantized even around E±_ = 0, i.e., at such large B (and therefore Aj, = 2), the disorder is not able to 
break the conductance quantization. As to the v = 4 state, which describes a large gap between the LLL and the 
n = 2 LL, its conductance is obviously quantized for all E±. 

B. Energy gaps and conductance: odd v 

The quantization of the conductance of the v = 1 state is broken at the same finite values of E± as those in the 
v = one^ It is also broken around E± = 0. As to the v — 3 state, the situation is similar to that in the v = 2 state: 
the quantization of its conductance is broken only around Ej_ =0. It is noticeable that in order to see these effects 
in the odd v states, one can use larger values of magnetic field (B = 3.43T in Ref. 0). This suggests that the odd v 
states arc more sensitive to disorder than those with even v. 

Let us discuss how these experimental data are reflected in the present analysis. The results for the v = 3 state, 
seems to be clear: as in the v = 2 one, in this case there is one phase describing the ground state and its gap has a 
minimum only at E± = 0. The case of the v = 1 state is more subtle. As for v = 0, there are two phases describing 
the ground state (the PSP and PLP ones). Moreover, the critical line dividing these phases exactly coincides with 
that between the SP and LP ones at v = 0. This is in accordance with the experimental data. However, unlike the 
v = state, there is no minimum in the gap at E" (see Fig. [3]). Instead, the gap is a smooth monotonically increasing 
function of E± and has a minimum only at E± = 0. This seems to suggest that, in disagreement with experiment, 
the quantization of the conductance at v = 1 should be broken only at E± = 0, despite the presence of the phase 
transition at E±_ = Ef ~ 5.04mV/nm. 

How can one solve this puzzle? A possible explanation might be as follows. Such a smooth behavior of the v = 1 
gap occurs due to delicate cancelations of different terms in calculating the gap in the PSP and PLP phases. This 
cancelation could be just an artifact of the approximation used in the present analysis. On the other hand, it may also 
suggest that although beyond this approximation a dynamics responsible for the conductance quantization breakdown 
at E±_ = Ef at v = 1 docs exist, it is weaker than that at v = 0, even in perfectly clean samples. This point deserves 
further study. 

VI. CONCLUSION 

In this paper, the analysis of Refs. Il6lli~7l was extended beyond the neutral point with the filling factor v = to 
describe the doped QH states with filling factors v = ±1, ±2, ±3, and v = ±4. It is noticeable that such a relatively 
simple model as the present one reproduces the main characteristics of experimentally observed broken-symmetry LLL 
states in a magnetic field. One important ingredient in the model is a very strong screening of the Coulomb interactions 
described by the polarization function. Such a strong screening radically changes the form of the interaction and is 
responsible for a linear scaling of dynamical gaps with a magnetic field in bilaycr graphene in contrast to monolayer 
graphene where a y/~B scaling takes place. The behavior of gaps as functions of magnetic field at zero electric field 
for different filling factors, FigO is found to be in a good agreement with experimental data£ The amplitudes of the 
gaps at the v = ±1,3 and v = ±2 plateaus are significantly smaller than the amplitude of the v = gap, due to 
the separate filling of the n — and n = 1 orbital Landau levels and the negative contribution of the Hartrce term, 
respectively. 

At nonzero electric field, we found a critical line in the plane of electric field, E±, and magnetic field, B, that 
separates both the spin and layer polarized phases at v = and the partially spin and layer polarized phases at 
v = ±1. On the other hand, there are unique phases for v = ±2 and v = ±3 QH states. The phase transition point 
moves out to larger electric fields as the magnetic field is increased, which implies that the ferromagnetic phase is 
stabilized by a magnetic field and is destabilized by an electric field. 

By studying the evolution of the gaps with an external electric field, we revealed a strong correlation between this 
evolution and the behavior of the conductance in the experiment^ the values of the electric field at the minima of the 
gaps correspond to those ones where the conductance is not quantized. The only exception is the v = ±1 QH state, 
where the gap has a minimum only at E± = and not at finite values of the electric field (possible reasons for such 
a disagreement were considered in Sec. IVBp . 
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Finally, we would like to note that as has been recently shown in Ref. l39l . nontrivial order parameters can also have 
important effects on the minimal dc conductivity in bilayer graphene. 

Note added: After this work was submitted for publication, a new experimental paper, Ref. Eol . has appeared, in 
which the phase diagram of the v = QH state in bilayer graphene is studied. In agreement with the theoretical 
predictions the critical line observed in that work separates the spin-polarized and layer-polarized states and has 
the form similar to the theoretical one, except the region with very low values of B, where the LLL approximation is 
apparently unreliable. In another recent experimental paper jil a phase transition to the nematic state was observed 
in bilayer graphene without magnetic field. The nematic state has been predicted in theoretical works42 It breaks 
the rotational symmetry and keeps quasiparticles gapless. It would be interesting to study a competition between the 
gapped state and the nematic state due to the Coulomb interaction and other interactions in bilayer graphene in a 
magnetic field. 
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Appendix A: Bare splitting of the Landau levels with n = and n = 1 

In this Appendix, we comment on the splitting of the Landau levels with orbital indices and 1 due to the bias 
electric field which leads to the sombrero shape of the upper band in the Hamiltonian of four-band model at zero 
magnetic ficldi^^ Under the reduction to the two-band model, the influence of the high-energy band is taken into 
account by the last term in the Hamiltonian ([2|). The bare splitting of the n = and n = 1 LLs is given by the 
formula^ 6 = edE±Hui c /-fi = A 5, where the dimensionless 5 is 6 = 2hw c /^i rj 0.011S[T]. 

This term leads to a modification of the gap equations (|22[) and (f2"3"| for the LLL with n = and n = 1: 

A -> (1 - nS)A . (Al) 
The free energy density (|27|) of the system will also receive a correction: 

£ -> £ = -i EE E [ E &»> + ^ + sZ - f^oC 1 - BgaC^en-) ■ (A2) 

£=± s=± n=0,l 

Note that the values of the magnetic field are less than 10T in the experiment with suspended bilayer graphene^ the 
results of our analysis are compared with. Therefore in this case the value of 5 = 0.011_B[T] is less than 0.1, and all 
these corrections are small. On the other hand, they could become relevant for stronger magnetic fields. 

By combining analytical and numerical methods, the modified gap equations and free energy density were analyzed. 
Now we list the main results that were obtained in this analysis. 

(1) For the v = state in Sec. IIV Al the transformation (|A2[) corresponds to E± in the free energy density of the 
LP phase effectively transforming as E± — >• (1 — S/2)E±. Because of that, the critical value E™ separating the SP 
and LP phases changes: 

The gaps for both the SP and LP phases are modified through the transformation E± (1 — S)E±: 

A^ - ^(h + h) + 2[Z (1 - S)eE ± d/2], A^ = - J(/ 2 + h) -2[Z+^-(l- 5)eE x d/2}. (A4) 
As a result, at the critical point, the two gaps will not coincide anymore. Instead, there is a jump in the gap there: 

A^ (Ao = A") - A^ (Ao = A"') = ^tj- 2 (z + ij?)- (A5) 
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(2) For the v = 1 state in Sec. IIVB1 the transformation (|A2[) corresponds to E± effectively transforming as E± —> 
(1 — S/2)E± in the free energy density in both the PSP and PLP phases. As a result, the critical value E" is modified 
in the same way as for the v = state and still coincides with the critical value in the latter. The gaps for the PSP 
and PLP solutions are modified as: 

Al^ = ^_(h+h-2I 2 )+5eE ± d/2, = -^(h + h - 2I 2 ) - 5eE±d/2. (A6) 

Therefore there is also a jump in the gap at the critical point: 

A^(A = A?) - A^(A„ = A"') = (z+tjs). (A7) 

It is noticeable that this jump coincides with that in the v = state (see Eq. (|A5[1 ). 

(3) For the v = 2 state, E± in the free energy density also effectively transforms as E± — > (1 — S/2)E±. There is a 
crossing of the energy levels £1,1,+ and E-i^i- in this state. As was shown in Sec. IIVC1 the value of Ef oas coincides 
with that of E" in the v = and v = 1 states, if the 5 correction is ignored. However, when this correction is 
included, these two values become different. Instead, at the crossing point, 

E" oss = -T, 77 \Z + • (A8) 

- 1 ed(l-5) \ kPJ v ' 

The form of the gap for E± < E" oss changes as 

ft 2 f 2e 2 dm\ 
A^ <s cro SS = (1 - S)eE±d + —p f I 2 + h - —fir ) ■ ( A9 ) 

The form of the gap for E± > E± oss remains unchanged. As a result, the kink singularity in the gap at the crossing 
point also does not change. 

(4) For the v = 3 state in Sec. IIVD1 the E± in the free energy density effectively transforms again as E± — > (l—5/2)E±. 
The gap in this phase changes to 

h 2 

A„ =3 = tt^2 ( J i + J 3 - 2/ 2) - SeE ± d/2. (A10) 
2ml z 
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